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1.  Introduction 

In  considering  transport  phenomena  in  biological  systems,  usually  the  purpose 
is  to  gain  some  information  about  the  structure  of  the  system  from  the  analysis 
of  flux  and  concentration  data.  In  a  typical  experiment  radioactively  labeled 
molecules  of  some  kind  are  introduced  into  an  animal ;  subsequently  radioactivity 
of  various  components  of  the  system  is  measured,  and  from  this  information  an 
attempt  is  made  to  characterize  the  transfer  of  molecules  from  one  part  of  the 
system  to  another  and  the  chemical  conversion  of  one  species  to  another.  Most 
commonly  this  characterization  has  been  in  terms  of  a  number  of  compartments 
or  “pools”  and  “turnover”  coefficients  between  the  different  pools.  A  recent 
review  article  by  Robertson  [1]  and  a  forthcoming  book  by  Sheppard  [2]  survey 
the  literature  in  this  area.  Mathematically  [3]  this  approach  can  be  shown  to 
represent  an  approximation  to  the  basic  transport  equation 

(1.1)  -div  J,(r,  t )  +  st(i,  t)  = 

where  J*(r,  t)  is  the  total  vector  flux  of  particles  of  the  kth  species  at  position  r 
and  time  t,  sk( r,  t)  is  the  net  production  per  unit  volume  from  chemical  reactions, 
and  Cfc(r,  t)  is  the  concentration.  Here  and  subsequently  we  always  refer  to  the 
labeled  particles  unless  specifically  stated  otherwise.  Physically  this  approach  can 
be  justified  by  the  existence  of  various  more  or  less  discrete  anatomical  and 
physiological  compartments  and  pools  in  biological  systems.  It  also  has  the 
practical  justification  that  frequently  flux  data  take  the  form  of  the  sum  of 
several  exponential  terms,  which  is  the  form  of  solution  obtained  for  the 
compartmental  model. 

The  main  difficulty  with  this  approach  is  that  it  quite  clearly  is  not  a  good 
approximation  for  certain  problems,  for  example  problems  in  which  transport 
via  the  circulatory  system  is  important.  Such  systems  can  of  course  be  described 
by  the  partial  differential  equation  (1.1),  but  this  is  essentially  vacuous  because 
the  equation  can  rarely  be  solved  for  biological  geometries.  An  attempt  has  been 
made  therefore  to  find  mathematical  descriptions  more  general  than  the  compart- 
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mental  approximation  and  yet  not  completely  useless  from  a  manipulative  point 
of  view.  One  such  is  to  describe  the  system  by  an  integral  equation  or  system  of 
integral  equations  in  which  the  kernel  is  regarded  as  characteristic  of  the  system. 
This  approach  has  been  applied  to  a  variety  of  problems.  (See  [3]  for  references.) 
Historically  the  method  was  applied  by  Volterra  to  a  variety  of  hereditary 
phenomena  including  biological  problems  [4],  [5].  In  earlier  papers  [3],  [6]  we 
have  discussed  in  detail  the  formulation  of  these  equations  for  certain  types  of 
biological  systems.  In  this  paper  we  will  consider  the  formulation  from  a  more 
general  point  of  view  and  will  discuss  the  problem  of  determining  the  kernel 
from  experimental  data.  To  illustrate  the  theory  we  will  briefly  discuss  a  partic¬ 
ular  problem  in  lipid  metabolism. 

2.  Formulation  of  equations 

We  will  assume  that  the  primary  experimental  variables  in  a  system  are  a 
set  of  fluxes,  finite  in  number,  denoted  by  71(f),  72(0 ,  •  •  •  ,  7 n(f).  For  the  moment 
we  will  consider  these  to  be  isotopically  labeled  molecular  fluxes,  although  the 
theory  we  develop  is  obviously  not  so  restricted.  We  make  the  fundamental 
assumption  that  each  flux  in  this  set  at  time  f  is  given  by  a  linear  functional  on 
all  previous  flux  in  the  set  plus  an  additive  term  accounting  for  material  being 
newly  introduced  into  the  system.  Specifically,  we  assume 

(2.1)  7.(0  =  E  P  co)7;(w)  do}  +  m<(f). 

i  J-°° 

If  we  give  the  fluxes  the  vector  representation 

(2.2)  rT  =  [71(0, 72(f),  •  •  • ,  7.(0], 

and  the  additive  term  the  representation, 

(2.3)  MT  =  [mi(t),  m2(t),  •  •  ,  mn{t)], 
the  system  (2.1)  can  be  written  as  the  operator  equation 

(2.4)  r  =  W  r  +  M. 

In  (2.4)  we  define  WJT  to  be  a  column  vector  with  fth  entry 

(2.5)  wr«(0  =  £  [*  Wij{t ,  w)7 fa)  do}. 

Equation  (2.4)  has  the  formal  solution 

(2.6)  r  =  (I  -  W)_1M 
and 

(2.7)  r  =  (I  +  W  +  W2  +  •••  +  W"  +  •  •  -)M. 

Either  (2.6)  or  (2.7)  gives  a  formal  solution  as  can  be  seen  by  direct  substi¬ 
tution  in  (2.4).  However,  the  series  in  (2.7)  must  be  computable  and  convergent 
to  have  any  real  meaning.  This  problem  will  be  discussed  below.  One  can  regard 
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equations  (2.6)  and  (2.7)  from  two  points  of  view,  depending  on  the  situation. 
One  is  that  the  operator  W  is  known  or  assumed  and  that  the  problem  is  to 
compute  the  response  of  the  system  to  an  arbitrary  input  M.  This,  for  example, 
is  the  problem  the  engineer  faces  in  computing  the  response  of  an  electrical 
network  composed  of  known  circuit  elements.  This  is  primarily  the  point  of 
view  taken  in  earlier  papers  [3],  [6].  The  other  point  of  view  is  to  regard  W  as 
unknown,  and  from  measurement  of  the  fluxes  after  known  input  to  determine  W. 
As  can  be  seen  from  (2.6),  what  can  actually  be  determined  experimentally  is 
(I  —  W)-1,  from  which  (I  —  W)  and  hence  W  can  be  found,  at  least  in  principle, 
by  some  appropriate  inversion  process.  It  is  this  problem  we  will  now  consider. 
Before  going  on  to  it,  however,  some  additional  remarks  can  be  made  about  the 
above  equations.  The  first  is  that  in  equation  (2.1)  the  summation  can  be  replaced 
by  integration  over  one  or  more  parameters,  say  over  space  or  momentum.  In 
actuality  fluxes  usually  are  more  or  less  continuously  distributed  in  space  and 
momentum,  but  in  biological  systems  this  problem  is  too  difficult  to  solve;  hence 
it  must  be  replaced  by  some  approximation  in  a  finite  number  of  fluxes.  Also, 
fluxes  integrated  over  some  range  of  space  and  momentum  usually  constitute 
the  actual  experimental  data.  The  second  remark  is  that  if  the  continuous  dis¬ 
tributions  in  time  are  replaced  by  discrete  terms,  the  vector  notation  and  the 
operator  formalism  remain  unchanged,  although  the  operator  instead  of  being 
of  the  integral  type  is  an  ordinary  matrix.  This  fact  is  highly  useful  in  numerical 
computations,  particularly  machine  computations,  and  also  frequently  is  of 
considerable  heuristic  value  in  more  general  considerations.  A  final  comment  is 
that  the  integrations  in  (2.1)  are  ordinarily  taken  from  some  initial  zero  time, 
marking  the  original  perturbation  of  the  system  by  the  introduction  of  labeled 
material. 


3.  Determination  of  the  operator  W  from  experimental  data 

3.1.  Finite  matric  representations.  We  will  first  consider  the  case  in  which  the 
continuous  time  distributuion  is  approximated  by  a  finite  number  of  discrete 
terms.  Specifically  we  will  suppose  that  in  (2.1)  each  y *(<)  is  replaced  by  a  finite 
number  of  discrete  numbers,  each  number  representing  the  total  flux  during  a 
time  interval  A U,  the  time  interval  of  integration  being  divided  into  m  such 
periods,  it  being  assumed  there  is  some  finite  lower  limit,  say  zero,  to  the  inte¬ 
gration.  Each  y i(t)  is  then  represented  by  a  column  vector  with  m  entries,  whose 
transpose  is 

(3.11)  7i  =  [7.1,  7.2,  •  *  •  ,  Tim]- 

Equation  (2.1)  then  becomes  the  linear  system 
(3.1.2)  7»i  =  2Z  Wijiky jk  H-  mu, 

j  k 

where  w,jik  is  the  probability  that  a  particle  which  contributes  to  the  jth  flux 
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in  the  kth  time  interval  contributes  to  the  fth  flux  in  the  lih  time  interval. 
Equation  (3.1.2)  can  be  written 

(3.1.3)  y .■  =  L  w ti7i  +  m„ 

j 

where  each  w,_,  is  an  m  by  m  matrix.  If  as  in  (2.2)  we  define 

(3.1.4)  rT  =  [71,  72,  •  •  *  ,  Yn]  =  [711,  •  ,  7lrn,  721,  •  •  •  ,  7 nm]  , 

then  r  is  a  column  matrix  with  nm  entries,  and  the  operator  equation  (2.4) 
retains  its  form 

(3.1.5)  r  =  wr  +  M, 

but  the  operator  W  is  now  a  matrix  built  up  from  the  Wq  matrices,  thus  with 
(nm)2  scalar  entries.  As  above  we  can  write  the  solution  of  (3.1.5) 

(3.1.6)  r  =  (I  -  W)-1M. 

Let  us  now  consider  the  problem  of  determining  (I  —  W)-1  from  experimental 
data.  Mathematically  it  is  obvious  that  (I  —  W)_I  is  completely  determined  by 
its  transformation  of  any  linearly  independent  set  of  inputs  nm  in  number,  that 
is,  any  basis  of  the  r  space.  Clearly  for  experimental  simplicity  the  orthonormal 
basis  (1,  0,  •  •  •  ,  0)  (0,  1,  0,  •  •  •  ,  0),  •  •  •  ,  (0,  •  ■  •  ,  0,  1)  is  the  logical  choice. 
Each  r  obtained  by  operating  on  this  basis  is  then  the  corresponding  column  of 
the  matrix  (I  —  W)-1.  That  is,  if  we  have 

(3.1.7)  Mf=  (1,0,  ...  ,0), 

(3.1.8)  A  =  (I  —  W)-iMlf 

then  the  entries  in  A  are  those  in  the  first  column  of  (I  —  W)-1  and  so  on.  Once 
(I  —  W)-1  is  known  it  can  be  inverted  to  give  (I  —  W)  and  hence  W.  We  assume 
it  is  nonsingular.  This  is  virtually  assured  by  the  physical  situation,  because 
if  not,  then  two  different  inputs  would  give  the  same  output,  or  to  put  it  dif¬ 
ferently  some  nonzero  input  would  give  zero  output.  Experimentally  the  vector 
(1,0,  •  •  •  ,0)  corresponds  to  the  introduction  of  unit  quantity  of  material  into 
the  system  with  momentum  and  spatial  distribution  corresponding  to  71  (t) 
during  the  first  time  period;  and  the  input  with  the  sth  entry  unity,  others 
zero,  and  s  —  (j  —  l)m  -f-  k,  corresponds  to  the  introduction  of  unit  material 
into  the  jth  flux  during  the  kth  time  period. 

3.2.  Volterra  systems.  Instead  of  the  above  finite  approximation  let  us  return 
to  the  system  of  integral  equations  (2.1).  The  operator  W  now  becomes  a  matrix 
composed  of  the  integral  operators  of  (2.1).  Thus,  by  definition,  we  have 

(3.2.1)  (IM)r  =  [mi(/),  m2(0,  •  •  •  ,  w„(0L 

(3.2.2)  (WM)r  =  [WMi(<),  WM2(0,  ,  WM,(<)], 

where 

(3.2.3)  WMj(f)  =  £  [‘  Wij(t,  co)m/co)  du. 

j  jo 
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As  indicated  by  (3.2.3)  WM  is  a  single  column  matrix  with  n  entries,  each  a 
function  of  time.  Successive  iterates,  which  appear  in  the  expansion  of 
(I  -  W^M,  are  defined  in  the  usual  way,  that  is, 

(3.2.4)  W2M  =  W(WM), 

(3.2.5)  W*M  =  W(WS_1M). 

Denoting  the  transpose  of  WS_1M  by 

(3.2.6)  (W^M)7  =  [W^M i(«),  •  •  •  ,  W-1Mn(0], 

we  obtain  for  the  ith  entry  in  W*M  from  the  operator  definitions 

(3.2.7)  W*M i  =  Z  [‘  co)Ws-1Mi(o))  doo. 

3 

The  resolvent  transformation  which  we  will  denote  by  K  is  then  defined  by 

(3.2.8)  (I  +K)(I  -  W)  =  I. 

Formally  we  clearly  have  for  K  the  expansion 

(3.2.9)  K  =  W  +  W2  +  •  •  •  +  Wn  +  •  •  •  , 
with  the  corresponding  kernel  matrix 

(3.2.10)  K (t,  w)  =  W (t,  «)  +  W *{t,  »)  +  •••+  W n(t,  «)+•••. 

The  entries  in  W (t,  go)  are  the  Wij(t,  go)  and  the  entries  in  W n(t,  go)  are  given  by 

(3.2.11)  w?j(t,  w)  =  L  T  wtk(t,  T)wntr\r,  go)  dr . 

k  Jo 

The  convergence  of  the  series 

(3.2.12)  KM  =  [W  +  W2  +  •  •  •  +  W»  +  •  •  -]M 

is  assured  by  the  fact  that  the  entries  in  W(£,  go)  are  Volterra  kernels,  that  is, 
W (t,  go)  =  0  if  go  >  t.  The  argument  remains  essentially  the  same  for  the  matrix 
as  for  a  single  function  (see  [7],  p.  147).  Briefly,  if  TFmax  is  the  upper  bound  of  the 

entries  in  W (t,  go)  and  A  is  the  upper  bound  of  the  integrals  th^go)  doo,  it  is 

easily  shown  that  the  entries  in  WM  are  maximized  by  wTFmaxA  and  the  entries 
in  W*M  by  (nTFmax)*A  £s-1/(s  —  1) !.  It  follows  that  the  series  for  KM  is  uniformly 
convergent  for  any  finite  upper  bound  of  t. 

Let  us  now  consider  the  problem  of  determining  the  operator  K  from  experi¬ 
mental  data.  We  have  the  general  operator  equation 

(3.2.13)  T  -  M  =  KM. 

We  will  suppose  that  at  some  time  go  unit  quantity  of  trace  material  is  suddenly 
(that  is,  as  rapidly  as  experimentally  possible)  introduced  into  the  jth.  flux,  all 
other  inputs  remaining  zero  throughout  the  experiment.  Mathematically  this 
corresponds  approximately  to  the  condition 

(3.2.14)  m3{t)  =  8(t  -  go). 

If  subsequently  the  fluxes  are  all  measured,  we  will  determine  the  entries  in 
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JT  —  M  for  this  injection  function.  The  fth  entry  in  the  column  matrix  JT  —  M 
will  then  determine  an  entry  co)  in  the  K (£,  co)  matrix.  In  other  words, 
with  this  isolated  experiment  we  can  determine  one  column  of  the  K (t,  co)  matrix 
for  one  value  of  co  and  all  subsequent  values  of  t.  This  is  essentially  all  the  infor¬ 
mation  which  can  be  obtained  from  a  single  experiment.  In  order  to  gain  addi¬ 
tional  information  further  separate  experiments  must  be  performed.  This 
obviously  requires  the  assumption  that  the  characteristic  matrix  for  the  system 
W (t,  co)  remains  the  same  in  the  separate  experiments.  In  other  words,  it  must 
be  possible  to  repeat  the  experiment  either  with  the  same  system  or  with  another 
nearly  identical  system.  Let  us  assume  that  this  can  be  done.  If  the  above  ex¬ 
periment  is  repeated,  with  the  injection  function  for  each  of  the  n  fluxes  being 
8(t  —  a>i)  sequentially,  each  experiment  will  determine  a  column  of  K (t,  coi),  and 
the  set  will  determine  the  matrix  K(t,  coi).  If  these  experiments  are  repeated  for 
a  set  of  values  of  co,  say  0,  coi,  •  •  •  ,  con,  the  corresponding  set  of  matrices  K (t,  0), 
K(£,  ooi),  •  •  •  ,  K {t,  co„)  will  be  obtained.  From  this  sequence  a  continuous  approxi¬ 
mation  of  K (t,  co)  which  will  be  designated  by  K'(f,  co)  can  be  constructed  by  some 
suitable  interpolation  scheme.  Then  from  K'(£,  co)  the  characteristic  matrix  of 
the  system  W(t,  co)  can  be  computed  approximately.  For  some  purposes  this 
may  be  superfluous,  the  entries  in  K '{t,  co)  giving  the  desired  information,  but 
usually  the  inversion  will  give  additional  information.  Basically,  the  inversion 
utilizes  the  general  operator  equation, 

(3.2.15)  (I  +  K)(I  -  W)  =  I, 
from  which  we  obtain  the  reciprocal  expansion 

(3.2.16)  W  =  K  -  K2  +K3  -  •••  -  (  —  l)nKn 

whose  convergence  is  assured  by  the  fact  that  K(£,  co)  is  again  a  Volterra  kernel. 
If  we  now  denote  K  —  K'  by  €,  substitution  gives 

(3.2.17)  W  =  (K'  +  e)  -  (K'  +  €)2  +  •  •  •  . 

Since  the  series  expansion  of  W  given  by  (3.2.16)  yields  a  uniformly  and  abso¬ 
lutely  convergent  series  when  operating  on  any  vector  for  which  the  integrals, 

/o'  dco,  are  bounded,  equation  (3.2.17)  can  be  rearranged  to  give 

(3.2.18)  W  =  K'  —  K/2  +  •  •  •  -  (-  l)"K'n  -  •  •  •  +  (terms  in  «). 

Here,  we  will  pass  over  a  discussion  of  the  terms  in  e,  that  is,  the  error  terms,  and 
simply  assume  they  are  small.  We  then  have 

(3.2.19)  W (t,  co)  ^  K '(*,  co)  -  K'2(J,  «)+•••-(-  l)nK'"(f,  co)  -  •  •  •  , 

where  the  iterations  have  the  same  definition  as  above. 

3.3.  Convolution  kernels.  It  can  be  seen  that  an  enormous  amount  of  ex¬ 
perimental  data  is  necessary  to  determine  a  general  Volterra  kernel  of  the  type 
W (t,  co).  Suppose  for  example,  that  one  has  a  modest  4  by  4  matrix,  and  takes  six 
time  values  co0,  coi,  co2,  co3,  co4,  co6.  This  requires  24  separate  experiments  for  one  com¬ 
plete  set  of  values.  To  test  for  significant  statistical  variation  in  the  entries,  sev- 
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eral  complete  sets  would  be  required.  In  systems  in  which  the  kernel  is  of  the 
convolution  type,  K (t  —  ui),  the  amount  of  experimental  data  necessary  for  an 
analysis  is  enormously  reduced,  because  each  column  need  be  determined  for  only 
one  value  of  c*>.  Thus,  for  a  4  by  4  matrix,  four  separate  experiments  would  give 
one  complete  set  of  data.  The  particular  time  «,  chosen  for  injection  clearly  is 
not  important.  Experimentally  systems  with  convolution  kernels  are  so-called 
“steady-state”  systems.  Once  K (t  —  oj)  is  determined  it  can  be  inverted  as  above 
If  K(f  —  co)  can  be  given  an  analytical  representation  the  mathematical  machin¬ 
ery  of  the  Laplace  transform  is  available  for  carrying  out  the  manipulations. 
A  detailed  analysis  of  systems  of  the  convolution  type  has  been  given  in  earlier 
papers  [3],  [6]  and  will  not  be  repeated  here.  It  is  worth  noting,  however,  that 
formally  the  analysis  of  these  systems  carries  over  to  systems  of  the  more  general 
type.  Hence,  the  operator  equations  developed  previously  for  these  systems 
remain  valid,  although  any  numerical  or  analytical  representation  becomes 
more  involved. 

3.4.  Systems  with  incomplete  information.  So  far  in  our  analysis  we  have  as¬ 
sumed  that  all  fluxes  are  available  for  injection  and  subsequent  measurement. 
This  may  not  be  the  case.  Let  us  suppose  that  s  are  available  for  injection,  t  for 
measurement.  This  will  permit  the  experimental  determination  of  st  entries  in 
the  JL(t,  w)  matrix.  Possibly  from  other  information  additional  entries  of  W (t,  co) 
or  of  K (t,  co)  are  known.  In  addition  some  of  the  unknown  entries  may  have  other 
restraints  on  them.  The  problem  is  to  utilize  all  the  available  information  to 
restrict  the  range  of  the  W (t,  co)  matrix.  The  prototype  of  this  problem  has  been 
considered  by  Berman  and  Schoenfeld  for  compartmental  systems.  In  our  dis¬ 
cussion  here,  we  will  do  little  more  than  state  the  general  problem. 

The  problem  appears  in  its  most  elementary  form  in  a  system  in  which  the 
operator  is  represented  by  a  finite  n  by  n  matrix.  As  before  we  have  the  general 
operator  equation 

(3.4.1)  (I  +  K)(I  -  W)  =  I, 

where  the  K  and  W  are  matrices.  If  the  multiplications  indicated  by  (3.4.1)  are 
carried  out,  we  will  obtain  n2  scalar  equations  in  the  2 n2  entries  in  (I  +  K)  and 
(I  —  W) .  Thus,  in  general,  n2  of  the  entries  can  be  expressed  as  functions  of  the 
other  n2.  If  these  are  all  known,  a  determinate  solution  is  obtained;  if  only  r  of 
them  are  known,  a  solution  with  n2  —  r  arbitrary  parameters  is  obtained.  Restric¬ 
tions  on  the  values  which  can  be  taken  by  these  parameters  and  the  solutions 
may  yield  useful  information. 

Another  way  of  looking  at  the  problem  is  to  consider  all  transformations 
P(I  +  K)  with  corresponding  transformation  of  the  inverse  (I  —  W)!*-1,  which 
will  preserve  the  known  restraints  on  K  and  W.  All  such  mappings  will  then 
generate  the  space  within  which  admissible  solutions  must  lie.  This  is  the  ap¬ 
proach  taken  by  Berman  and  Schoenfeld  in  their  paper  on  compartmental 
systems  [8]. 

The  problem  can  quite  clearly  be  extended  to  integral  operators  of  the  Volterra 
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type,  or  for  that  matter  to  linear  operators  in  general.  As  yet  the  only  case  we 
have  investigated  at  all  systematically  is  that  in  which  the  integral  operators 
are  of  the  convolution  type.  Here,  instead  of  (3.4.1)  we  have 

(3.4.2)  [I  +  K(p)][I  -  W (p)]  =  I, 

where  K(p)  and  W(p)  are  matrices  composed  of  the  Laplace  transforms  of  the 
entries  in  K(f)  and  W (t).  The  remarks  made  above  apply  except  that  the  numeri¬ 
cal  entries  in  the  K  and  W  matrices  are  replaced  by  the  corresponding  functions 
in  the  transform  parameter  p.  Hence,  we  will  have  n2  relations,  involving  2n2 
functions.  From  these  equations  n2  of  the  functions  can  be  expressed  (at  least 
theoretically)  in  terms  of  the  remaining  n2.  If  r  of  these  are  known,  we  will 
obtain  a  solution  in  terms  of  n2  —  r  arbitrary  functions  of  p.  As  before  rather 
general  restraints  may  restrict  the  range  of  possible  solutions  and  lead  to  useful 
information. 

4.  Application  to  a  kinetic  problem 

4.1.  Analysis  of  fatty  acid  transport.  For  obvious  reasons,  steady-state  sys¬ 
tems,  that  is,  systems  with  a  kernel  of  the  convolution  type,  have  received  more 
practical  attention,  and  the  above  theory  will  be  illustrated  with  a  rather  simple 
application  of  this  type.  Nevertheless,  conceptually  the  analysis  of  a  much  more 
complicated  system  would  be  much  the  same. 

Our  illustrative  application  arose  in  the  experimental  studies  of  Dr.  D.  S. 
Fredrickson  and  colleagues  on  lipid  transport  and  metabolism.  Basically  the 
question  is  whether  unesterified  fatty  acid  (denoted  by  UFA)  and  fatty  acid 
derived  from  the  splitting  of  triglyceride  (TGFA)  follow  the  same  metabolic 
pathway — more  exactly  whether  the  fatty  acid  derived  from  splitting  triglyc¬ 
eride  all  or  nearly  all  passes  through  the  plasma  pool  of  unesterified  fatty  acid 
before  it  is  metabolized.  The  experimental  details  of  the  work  have  been  de¬ 
scribed  previously  [9]  and  here  we  give  only  an  outline  of  the  mathematical 
analysis,  the  details  of  which  will  be  published  elsewhere. 

Two  sets  of  experiments  were  performed.  In  one  set,  C14  labeled  UFA  was 
injected  directly  into  the  blood  plasma  of  dogs.  Subsequently  the  specific  activity 
of  UFA  in  the  plasma,  and  the  specific  activity  of  C14  labeled  carbon  dioxide 
in  the  respiratory  output  were  measured.  In  the  other  set  C14  labeled  TGFA  was 
injected  into  the  plasma  as  chylomicra.  Subsequently  the  specific  activities  of 
plasma  UFA  and  of  respired  CO2  were  measured.  In  terms  of  these  experiments 
the  problem  was  to  determine  whether  in  the  second  set,  the  flux  of  labeled 
UFA  out  of  the  plasma  was  sufficient  to  account  for  the  total  output  of  labeled 
C02  or  whether  a  significant  fraction  of  TGFA  was  metabolized  without  passing 
through  the  plasma  UFA  pool.  To  decide  this  question  it  is  clearly  necessary  to 
compute  the  output  of  labeled  CO2  from  the  labeled  UFA  which  passes  through 
the  plasma  pool.  This  can  be  done  utilizing  the  above  mathematical  theory  and 
information  from  the  first  set  of  experiments. 
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To  set  up  our  model  we  will  suppose  that  we  have  two  compartments,  plasma 
and  tissue.  The  plasma  is  considered  to  be  uniformly  mixed,  but  the  tissue  is 
not.  A  labeled  UFA  molecule  leaving  the  plasma  has  three  possible  fates: 
(1)  after  a  certain  interim  it  returns  to  the  plasma  as  an  UFA  molecule,  (2)  it  is 
metabolized  and  its  C14  appears  in  the  respired  CO2,  (3)  it  is  deposited  in  some 
tissue  depot  where  for  practical  purposes  it  remains  indefinitely  (of  the  order  of 
months  or  years).  It  must  be  recognized,  of  course,  that  en  route  to  these  fates 
the  UFA  molecules  may  follow  complicated  and  devious  pathways  through 
other  lipid  compartments  located  in  both  tissues  and  plasma.  Let  us  designate 
the  flux  of  UFA  molecules  leaving  the  plasma  by  71(f),  the  returning  flux  by 
72(f),  and  the  flux  of  labeled  C02  by  73(f).  We  assume  that  these  fluxes  are 
related  by  integral  operators  of  the  Yolterra  type,  and  the  experiments  were 
performed  so  that  the  dogs  were  in  as  near  a  physiological  steady  state  as 
possible.  Hence,  it  is  reasonable  to  assume  that  these  operators  are  of  the  con¬ 
volution  type.  Thus,  we  obtain  the  equations  for  the  system 

(4.1.1)  72(f)  =  Jo  w2i(t  —  00)71(01)  do), 

(4.1.2)  73  (f)  =  Jq  w3i(t  —  01)71(01)  do). 

In  the  first  set  of  experiments  the  fluxes,  71(f),  72(f),  and  73(f),  are  known  from 
the  experimental  data:  73(f)  is  measured  directly  and  71(f)  and  72(f)  are  deter¬ 
mined  from  measurements  of  plasma  UFA  activity,  denoted  by  Ci(f),  in  a  way 
we  will  now  briefly  outline.  If  V  is  the  plasma  volume,  we  have  for  the  total 
quantity  of  labeled  UFA  in  the  plasma 

(4.1.3)  gi(f)  =  ci(f)F. 

The  rate  of  change  of  the  total  quantity  of  labeled  UFA  in  the  plasma  pool 
obviously  equals  the  difference  between  the  ingoing  and  outgoing  fluxes,  or 

(4.1.4)  ^  =  7 ,«)  -  7.(1)- 

During  the  initial  phase  of  the  experiment,  UFA  radioactivity  in  the  plasma 
falls  very  rapidly  and  nearly  exponentially.  It  is  reasonable  to  assume  that  during 
this  initial  period  72(f)  =  0.  With  this  assumption  and  the  assumption  that  the 
plasma  is  well  mixed,  we  have 

(4.1.5)  7.(0)  =  kcMV  =  kqt( 0)  = 

where  k  is  a  turnover  coefficient,  which  can  be  computed  from  (4.1.5).  Subse¬ 
quently  we  assume  that  k  remains  constant,  that  is,  that  a  fraction  k  of  labeled 
UFA  molecules  present  in  the  plasma  UFA  pool  leave  it  per  unit  time,  so  that 

(4.1.6)  7i(f)  =  kqi(t)  =  kci(t)V, 

in  which  k,  ci(f)  and  V  are  known  directly  from  the  experimental  data  [V  = 
m/qi(0),  where  m  is  the  total  quantity  of  labeled  UFA  introduced  initially].  The 
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flux  of  UFA  entering  the  plasma,  72 (t),  can  then  be  computed  from  (4.1.4). 
Knowing  the  fluxes,  yi(t),  72 (t)  and  73 (t),  we  can  determine  Wn(t  —  w)  and 
Wn(t  —  w)  by  solving  equations  (4.1.1)  and  (4.1.2).  We  have  done  this  analyti¬ 
cally  by  approximating  7 x(<)  and  72 (t)  and  yz(t)  by  sums  of  exponentials  and 
using  Laplace  transforms.  We  have  also  carried  out  numerical  solutions  on  an 
IBM  650.  Both  w2i(t  —  w)  and  w3i(t  —  u)  have  physiological  significance,  but 
it  is  the  latter  we  want  for  our  computation  of  the  labeled  CO2  due  to  plasma 
UFA  in  the  second  set  of  experiments. 

In  the  second  set  of  experiments,  the  chylomicra  are  removed  from  the  plasma 
and  enter  the  tissue.  Here,  they  are  at  least  partially  split  into  fatty  acids  and 
glycerol;  at  least  a  fraction  of  the  fatty  acid  returns  to  the  plasma  and  enters 
the  plasma  UFA  pool.  Subsequently  this  labeled  UFA  leaves  the  plasma  with 
the  possible  fates  given  above.  If  we  assume  that  the  turnover  constant  k  is  the 
same  as  in  the  first  set  of  experiments,  a  condition  which  is  reasonably  assured  by 
the  equivalent  total  UFA  concentration  in  the  two  sets  of  experiments  [10],  71  (t) 
can  be  computed  from  the  plasma  level  of  labeled  UFA  as  in  the  earlier  experi¬ 
ments,  that  is,  by  means  of  (4.1.6) ;  and  if  the  transport  function  w3i(t  —  a>)  is  as¬ 
sumed  to  remain  the  same,  the  production  of  labeled  CO2  from  71  (t)  can  be 
computed  by  equation  (4.1.2).  The  total  production  of  labeled  CO2  is  measured 
directly.  The  excess  of  the  total  production  of  C02  over  that  computed  to  arise 
from  71(f),  gives  the  production  from  TGFA  whose  fatty  acids  do  not  cycle 
through  the  plasma  pool  of  UFA.  We  have  not  yet  completed  our  analysis  of 
these  experiments,  but  it  turns  out  that  a  considerable  fraction  of  TGFA  is 
oxidized  without  passing  through  the  plasma  UFA  pool. 

5.  Concluding  remarks 

In  essence  in  the  theory  presented  in  this  paper  we  represent  the  state  of  a 
biological  system  by  a  vector  in  an  abstract  space  and  consider  the  basic  char¬ 
acterization  of  the  system  to  be  a  linear  operator  in  the  space.  We  have  outlined 
how  the  structure  of  the  operator  can  be  obtained  from  experimental  data.  It 
is  clear  that  this  representation  is  primarily  a  data  reduction  scheme.  Until  the 
method  has  been  applied  to  a  larger  number  of  biological  problems,  it  is  difficult 
to  assess  its  ultimate  usefulness.  One  of  the  greatest  assets  of  this  approach  is 
that  it  provides  a  precise  mathematical  language  in  which  a  large  number  of 
systems  can  be  described  and  in  terms  of  which  experiments  can  be  formulated. 
This  we  have  tried  to  illustrate  with  the  example  from  fatty  acid  transport  and 
metabolism.  Another  great  advantage  of  the  method  is  that  both  stationary  and 
nonstationary  problems  can  be  handled  by  it,  utilizing  the  same  formal  language. 
The  method  does  not  yield  any  direct  information  about  the  microscopic  struc¬ 
ture  of  a  system,  but  any  microscopic  theory  must  yield  an  operator  which 
agrees  with  that  derived  directly  from  the  experimental  data,  in  much  the  same 
way  that  thermodynamic  coefficients  computed  by  statistical  theory  must  agree 
with  those  experimentally  measured. 
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It  is  apparent  that  we  have  raised  a  number  of  unanswered  problems  in  this 
paper.  There  are  statistical  sampling  problems  relative  to  the  construction  of 
the  operator  from  experimental  data.  We  have  just  touched  on  the  question  of 
the  analysis  of  systems  with  incomplete  information.  We  hope  that  these  and 
other  problems  will  be  solved  by  further  investigation  along  the  directions  we 
have  outlined. 

The  author  thanks  Mr.  Arnold  Jones  for  his  great  help  in  preparing  this  paper. 
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